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£Sj ■ Neutron stars are some of the densest manifestations of massive objects in the 



universe. They are ideal astrophysical laboratories for testing theories of dense 
matter physics and provide connections among nuclear physics, particle physics 
and astrophysics. Neutron stars may exhibit conditions and phenomena not ob- 
served elsewhere, such as hyperon-dominated matter, deconfined quark matter, 
^ ■ superfluidity and superconductivity with critical temperatures near 10 10 kelvin, 

^ ! opaqueness to neutrinos, and magnetic fields in excess of 10 13 Gauss. Here, we 

describe the formation, structure, internal composition and evolution of neutron 
O ■ stars. Observations that include studies of binary pulsars, thermal emission from 



Or 



isolated neutron stars, glitches from pulsars and quasi-periodic oscillations from 



Q ■ accreting neutron stars provide information about neutron star masses, radii, 

; 

c/2 ! temperatures, ages and internal compositions. 
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j ■ »-j ■ The term neutron star as generally used today refers to a star with a mass M on the 

order of 1.5 solar masses (M ), a radius R of ~ 12 km, and a central density n c as high as 5 
to 10 times the nuclear equilibrium density n ~ 0.16 fm -3 of neutrons and protons found in 



laboratory nuc 
universe 



ei. A neutron star is thus one of the densest forms of matter in the observable 
l], 0, • Although neutrons dominate the nucleonic component of neutron stars, 
some protons (and enough electrons and muons to neutralize the matter) exist. At supra- 

nn 

nuclear densities, exotica such as strangeness-bearing baryons |4],|5j], condensed mesons (pion 
or kaon) 0, Q, la], or even deconfined quarks P| may appear. Fermions, whether in the 
form of hadrons or deconfined quarks, are expected to also exhibit superfluidity and/or 
superconductivity. 

Neutron stars encompass "normal" stars, with hadronic matter exteriors in which the 
surface pressure and baryon density vanish (the interior may contain any or a combination of 
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exotic particles permitted by the physics of strong interactions), and "strange quark matter" 
(SQM) stars [lCj. An SQM star could have either a bare quark matter surface with vanishing 
pressure but a large, supra- nuclear density, or a thin layer of normal matter supported by 
Coulomb forces above the quark surface. The name SQM star originates from the conjecture 
that quark matter with up, down and strange quarks (the charm, bottom and top quarks are 
too massive to appear inside neutron stars) might have a greater binding energy per baryon 
at zero pressure than iron nuclei have. If true, such matter is the ultimate ground state 
of matter. Normal matter is then metastable, and compressed to sufficiently high density, 
would spontaneously convert to deconfined quark matter. Unlike normal stars, SQM stars 
are self-bound, not requiring gravity to hold them together. It is generally assumed that 
pulsars and other observed neutron stars are normal neutron stars. If SQM stars have a 
bare quark surface, calculations suggest suggest that photon emission from SQM stars occurs 



primarily in the energy range 30 keV < E < 500 keV 



ill- 



How Neutron Stars are Formed 



Neutron stars are created in the aftermath of the gravitational collapse of the core of a 
massive star (> 8 M Q ) at the end of its life, which triggers a Type II supernova explosion. 
Newly-born neutron stars or proto-neutron stars are rich in leptons, mostly e~ and v e (Fig. 
1). The detailed explosion mechanism of Type II supernovae is not understood [3], but it 
is probable that neutrinos play a crucial role. One of the most remarkable aspects is that 
neutrinos become temporarily trapped within the star during collapse. The typical neutrino- 
matter cross section is cr ~ 10~ 40 cm 2 , resulting in a mean free path A ~ (em)" 1 « 10 cm, 
where the baryon number density is n ~ 2 to 3 n . This length is much less than the 
proto-neutron star radius which exceeds 20 km. The gravitational binding energy released 
in the collapse of the progenitor star's white dwarf-like core to a neutron star is about 
3GM/5R 2 ~ 3 x 10 53 erg (G is the gravitational constant), which is about 10% of its total 
mass energy Mc 2 . The kinetic energy of the expanding remnant is on the order of 1 xlO 51 to 
2 xlO 51 erg, and the total energy radiated in photons is further reduced by a factor of 100. 
Nearly all the energy is carried off by neutrinos and antineutrinos of all flavors in roughly 
equal proportions. 

Core collapse halts when the star's interior density reaches no, which triggers the for- 
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mation of a shock wave at the core's outer edge. The shock wave propagates only about 
100 to 200 km before it stalls, having lost energy to neutrinos and from nuclear dissociation 
of the material it has plowed through [stage (I) in Fig. 1]. Apparently, neutrinos from the 
core, assisted perhaps by rotation, convection and magnetic fields, eventually resuscitate 
the shock, which within seconds accelerates outwards, expelling the massive stellar mantle. 
The proto-neutron star left behind rapidly shrinks because of pressure losses from neutrino 
emission in its periphery (stage II). The escape of neutrinos from the interior occurs on a 
diffusion time r ~ 3R 2 /Xc ~ 10 s. The neutrinos observed from Supernova (SN) 1987A 
in the Large Magellanic Cloud confirmed this time scale and the overall energy release of 



~ 3 x 10 5d ergs 



14, 
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16]. 



The loss of neutrinos (which forces electrons and protons to combine, making the matter 
more neutron-rich) initially warms the stellar interior. The core temperature more than 
doubles (stage III), reaching ~ 50 MeV (6 x 10 11 K). After 10 to 20 s, however, the steady 
emission of neutrinos begins to cool the interior. Because the cross section a oc A -1 scales 
as the square of the mean neutrino energy, the condition A > R is achieved in about 50 s. 
The star becomes transparent to neutrinos (stage IV), and its cooling rate accelerates. 

Neutron stars have both minimum and maximum mass limits. The maximum mass, 
which is of purely general relativistic origin, is unknown, but lies in the range of 1.44 to 3 
M . The upper oonnd follows frona ean S aHt y Q, that the speed of sonnd . dense .natter 
is less than the speed of light, whereas the lower bound is the largest accurately measured 
pulsar mass, 1.4408 ± 0.0003 M Q , in the binary pulsar PSR 1913+16 jisl ]. The minimum 
stable neutron star mass is about 0.1 M , although a more realistic minimum stems from a 
neutron star's origin in a supernova. Lepton-rich proto-neutron stars are unbound if their 
masses are less than about 1 M Q jl^ . 

The proto-neutron star, in some cases, might not survive its early evolution, collapsing 
instead into a black hole. This could occur in two different ways. First, proto-neutron 
stars accrete mass that has fallen through the shock. This accretion terminates when the 
shock lifts off, but not before the star's mass has exceeded its maximum mass. It would 
then cohapse and its nent.no signa. „on, d a bfflP t. y eease Q. If this does not occur, a 



2l| . A proto-neutron star's maximum mass 



second mode of black hole creation is possible 
is enhanced relative to a cold star by its extra leptons and thermal energy. Therefore, 
following accretion, the proto-neutron star could have a mass below its maximum mass, but 
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still greater than that of a cold star. If so, collapse to a black hole would occur on a diffusion 
time of 10 to 20 s, longer than in the first case. Perhaps such a scenario could explain the 
enigma of SN 1987A. The 10 s duration of the neutrino signal confirmed the birth and 
early survival of a proto-neutron star, yet there is no evidence that a neutron star exists 
in this supernova's remnant. The remnant's observed luminosity is fully accounted for by 
radioactivity in the ejected matter 2^ . meaning that any contribution from magnetic dipole 
radiation, expected from a rotating magnetized neutron star, is very small. Either there is 
presently no neutron star, or its spin rate or magnetic field are substantially smaller than 
those of typical pulsars. A delayed collapse scenario could account for these observations j^lj . 

Global Structure of Neutron Stars 

Global aspects of neutron stars, such as the mass-radius (M—R) relation, are determined 

relativity (GR), 



24] 



by the equations of hydrostatic equilibrium. For a spherical object in genera 
these are the so-called TOV (Tolman-Oppenheimer-Volkov) equations 0, 

dP G(m(r) + 4vrr 3 P/c 2 ) (p + P/c 2 ) 
dr r(r — 2Gm(r)/c 2 ) 

= iW, (i) 

where P and p are the pressure and mass-energy density, respectively, and m(r) is the 
gravitational mass enclosed within a radius r. Although a few exact solutions are known 
|25| . for a realistic P — p relation (equation of state, hereafter EOS) these equations must 
be numerically solved to obtain the M — R relation as shown in Fig. 2. The region in Fig. 
2 bounded by the Schwarzschild condition R < 2GM/c 2 is excluded by general relativity, 
and that bounded by R < 3GM/c 2 is excluded by causality j^J. Some normal neutron star 
cases, such as GS1, contain large amounts of exotica, any of which produces a large amount 
of softening and relatively small radii and maximum masses. For small masses, SQM stars 
are nearly incompressible (R oc M 1 / 3 ). 

For normal neutron stars the radius is relatively insensitive to the mass in the vicinity of 
1 to 1.5 M Q unless the maximum mass is relatively small. A simultaneous measurement of 
mass and radius of an intermediate mass star could help to discriminate among the families of 
possible EOS's. Perhaps two of the most important, but unknown, astrophysical quantities 
are the neutron star maximum mass and the radius of 1.4 M neutron stars. 



There are large variations in predicted radii and maximum masses (Fig. 2) because of 
the uncertainties in the EOS near and above no [27]. This seems paradoxical because the 
properties of matter inside laboratory nuclei are thought to be well understood. However, 
an important distinction between nuclear and neutron star matter is their relative proton 
fraction x. Nuclear matter has nearly equal numbers of neutrons and protons (x ~ 1/2), 
but neutron star matter has only a few percent protons. The energy can be described with 
a quadratic interpolation in the proton fraction x : 

E(n,x) = E{n,x= 1/2) + S v {n)(l - 2xf . (2) 

The symmetry energy function S v (n) is uncertain, although weak constraints exist from 
ground-state masses (binding energies) and giant dipole resonances of laboratory nuclei. 
The symmetry energy of nuclei is divided between bulk and surface contributions, which 
scale with nuclear mass number as A and A 2 / 3 , respectively, but the ranges of A 1 / 3 (up to 
6) and x in laboratory nuclei are too small to separate them. 

A consequence of this uncertainty is that different models predict up to a factor of 
6 variation in the pressure of neutron star matter near no, even though the pressure of 
symmetric matter is better known, being nearly zero at the same density. This pressure 
variation accounts for the nearly 50% variation in predictions of neutron star radii j^. 

A potential constraint on the EOS derives from the rotation of neutron stars. An 
absolute upper limit to the neutron star spin frequency is the mass-shedding limit, at which 
the velocity of the stellar surface equals that of an orbiting particle suspended just above 
the surface. For a rigid Newtonian sphere this frequency is the Keplerian rate 

VK = {2n)- l ^GM/R 3 = 1833 (M/M ) 1/2 (10 km/R) 3/2 Hz . (3) 

However, both deformation and GR effects are important. A similar expression, but with a 
coefficient of 1224 Hz and in which M and R refer to the mass and radius of the maximum 
mass, non-rotating, configuration, describes the maximum rotation rate possible for an EOS 
|2fiL I28L l2Ql| . We have found that Eq. (J3}, but with a coefficient of 1045 Hz, approximately 
describes the maximum rotation rate for a star of mass M (not close to the maximum mass) 
and non-rotating radius R independently of the EOS. The highest observed spin rate, 641 
Hz from pulsar PSR B1937+21 j^jj], implies a radius limit of 15.5 km for 1.4 M . 
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Internal Structure and Composition 



A neutron star has five major regions, the inner and outer cores, the crust, the envelope 
and the atmosphere (Fig. 3). The atmosphere and envelope contain a negligible amount of 
mass, but the atmosphere plays an important role in shaping the emergent photon spectrum, 
and the envelope crucially influences the transport and release of thermal energy from the 
star's surface. The crust, extending about 1 to 2 km below the surface, primarily contains 
nuclei. The dominant nuclei in the crust vary with density, and range from 56 Fe for matter 
with densities less than about 10 6 g cm -3 to nuclei with A ~ 200 but x ~ (0.1 to 0.2) near 
the core-crust interface at n m n /3. Such extremely neutron-rich nuclei are not observed 
in the laboratory, but rare-isotope accelerators jjil] hope to create some of them. 

Within the crust, at densities above the neutron drip density 4 x 10 11 g cm -3 where 
the neutron chemical potential (the energy required to remove a neutron from the filled 
sea of degenerate fermions) is zero, neutrons leak out of nuclei. At the highest densities 
in the crust, more of the matter resides in the neutron fluid than in nuclei. At the core- 
crust interface, nuclei are so closely packed that they are almost touching. At somewhat 
lower densities, the nuclear lattice can turn inside-out and form a lattice of voids, which 
is eventually s q uee Z ed out at densities nea r „ Q. If so, beginning at about 0, n „, tbe r e 
could be a continuous change of the dimensionality of matter from three-dimensional (3-D) 
nuclei (meatballs), to 2-D cylindrical nuclei (spaghetti), to 1-D slabs of nuclei interlaid with 
planar voids (lasagna), to 2-D cylindrical voids (ziti) , to 3-D voids (ravioli, or Swiss cheese 
in Fig. 3) before an eventual transition to uniform nucleonic matter (sauce). This series of 
transitions is known as the nuclear pasta. 

For temperatures less than ~ 0.1 MeV, the neutron fluid in the crust probably forms a x Sq 
superfluid |l|. Such a superfluid would alter the specific heat and the neutrino emissivities 
of the crust, thereby affecting how neutron stars cool. The superfluid would also form a 
reservoir of angular momentum that, being loosely coupled to the crust, could cause pulsar 
glitch phenomena jsjj. 

The core constitutes up to 99% of the mass of the star (Fig. 3). The outer core consists 
of a soup of nucleons, electrons and muons. The neutrons could form a 3 P2 superfluid and 
the protons a 1 So superconductor within the outer core. In the inner core exotic particles 
such as strangeness-bearing hyperons and/or Bose condensates (pions or kaons) may become 
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abundant. It is possible that a transition to a mixed phase of hadronic and deconfined quark 

timate ground state of matter, 
has yielded novel states of 



matter develops [34J , even if strange quark matter is not the u 
Delineating the phase structure of dense cold quark matter 



matter including color-superconducting phases with 



and without condensed mesons 13 



Neutron Star Cooling 

The interior of a proto-neutron star loses energy at a rapid rate by neutrino emission. 
Within 10 to 100 years, the thermal evolution time of the crust, heat transported by electron 
conduction into the interior, where it is radiated away by neutrinos, creates an isothermal 
structure (stage (V) in Fig. 1). The star continuously emits photons, dominantly in x-rays, 
with an effective temperature T e ff that tracks the interior temperature but that is smaller 
by a factor of ~ 100. The energy loss from photons is swamped by neutrino emission from 
the interior until the star becomes about 3 x 10 5 years old (stage VI). 

The overall time that a neutron star will remain visible to terrestrial observers is not yet 
known, but there are two possibilities: the standard and enhanced cooling scenarios. The 



dominant neutrino cooling reactions are of a general type, known as Urea processes |37j, in 
which thermally excited particles alternately undergo beta and inverse-beta decays. Each 
reaction produces a neutrino or antineutrino, and thermal energy is thus continuously lost. 
The most efficient Urea process is the direct Urea process involving nucleons: 

n — > p + e~ + v e , p — > n + e + + v e . (4) 

This process is only permitted if energy and momentum can be simultaneously conserved. 
This requires that the proton to neutron ratio exceeds 1/8, or the proton fraction x > 1/9, 
which is far above the value found in neutron star matter in the vicinity of no- In a mixture 
of neutrons, protons and electrons, the proton fraction x in beta equilibrium satisfies j^j 

x ~ 0.048 {S v (n)/S v (n )) s (n /n) (1 - 2a;) 3 , (5) 

where, typically, S v (n ) ~ 30 MeV. Because x generally increases with density, the direct 
Urea process might still occur above some density threshold. However, if the direct process 
is not possible, neutrino cooling must occur by the modified Urea process 

n + (n, p) — > p + (n, p) + e~ + v e , p + (n,p) — > n + (n,p) + e + + v e , (6) 



in which an additional nucleon (n,p) participates in order to conserve momentum. The 
modified Urea rate is reduced by a factor of (T//j, n ) 2 < 10 -4 to 10~ 5 compared to the 
direct Urea rate, and neutron star cooling is correspondingly slower. The standard cooling 
scenario assumes that direct Urea processes cannot occur, and predicts that neutron stars 
should remain observable by surface thermal emission for up to a few million years. 

The question of whether or not the direct Urea process occurs in neutron stars is of 
fundamental importance. The density dependence of the symmetry energy function S v 
determines the values of x and the threshold density at which the nucleonic direct Urea 
process occurs (Eq. ©). It also plays an essential role in determining the threshold densities 
of other particles, such as hyperons, pions, kaons or quarks, whose existences trigger other 
direct Urea processes 3^|. If a star's central density lies below the Urea threshold, enhanced 
cooling cannot occur. Again, the quantity S v (n) plays a crucial role for neutron stars, and 
its inherent uncertainty means that it is presently unknown if direct Urea processes can 
occur in neutron stars. 

There are two additional issues affecting cooling trajectories of neutron stars: superflu- 
idity 2i3, Q| and envelope composition . Superfluidity quenches cooling from the direct 
Urea process. However, an additional cooling source from the formation and breaking of 
nucleonic Cooper pairs increases the cooling rate from the modified Urea process Nev- 
ertheless, a clear distinction remains between enhanced and standard cooling trajectories. 

Envelope composition also plays a role in the inferred surface temperatures. Although 
it is commonly assumed that the envelope is dominated by iron-peak nuclei, this may not 
be the case. Light elements (H or He) have smaller photon opacities which enhance surface 
photon emission. Neutron stars appear warmer with light-element envelopes for their first 
100,000 years of cooling but eventually the situation reverses 

Observations and Inferred Stellar Properties 



Masses 



The most accurately measured neutron star masses are from timing observations of radio 

n 

binary pulsars |44|. These include pulsars orbiting another neutron star, a white dwarf or 
a main-sequence star. Ordinarily, observations of pulsars in binaries yield orbital sizes and 
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periods from Doppler shift phenomenon, from which the total mass of the binary can be 
deduced. But the compact nature of several binary pulsars permits detection of relativistic 
effects, such as Shapiro delay j3] or orbit shrinkage due to gravitational radiation reaction, 
which constrains the inclination angle and permits measurement of each mass in the binary. 
A sufficiently well-observed system, such as the binary pulsar PSR 1913+16 or the newly 
discovered pulsar binary PSR J0737-3039 4fJ, can have masses determined to impressive 
accuracy. Masses can also be estimated for neutron stars that are accreting matter from 
a stellar companion in so-called x-ray binaries, but the measurements have much larger 
relative errors Table |U Neutron stars in binaries with white dwarf companions have a 
broader range of masses than binary neutron stars and the wider mass range may signify 
a wider range of formation mechanisms. It has been suggested that a rather narrow set of 
evolutionary circumstances conspire to form double neutron star binaries j^. The largest 
apparent masses are in the systems 4U1700-37, which might in fact contain a black hole, 
not a neutron star, Vela X-l, and the pulsar J0751+1807, but all have large uncertainties. 
Raising the limit for the neutron star maximum mass could eliminate entire EOS families, 
especially those in which exotica appear and substantial softening begins around 2 to 3 no- 
This could be significant, because exotica generally reduce the maximum mass appreciably. 

Thermal Emission 

Most known neutron stars are observed as pulsars and have photon emissions from radio 
to x-ray wavelengths dominated by non-thermal emissions. It is believed that the bulk of 
the non-thermal emissions are generated in a neutron star's magnetosphere. Although such 
emissions can teach us about magnetospheric phenomena, they are difficult to utilize in 
constraining the star's global aspects, such as mass, radius and temperature, that have 
significant bearing on a star's interior structure, composition and evolution. About a dozen 
neutron stars with high thermal emissions, and with ages up to a million years, have been 



identified |43( and these stars are expected in the standard cooling scenario to have surface 
temperatures in the range of 3 x 10 5 to 10 6 K (Fig. 4), so the bulk of their emitted radiation 
should lie in the extreme ultraviolet or x-ray regions. 
The effective temperature T e ff )OQ is defined from 

= LjAitd 2 = a B T* fLoo (R^/df , (7) 
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where <7b is the Stefan-Boltzmann constant, d is the distance, and Foo and refer to the 
flux and luminosity observed at Earth. These latter quantities, and T e //,oo, are redshifted 
from the neutron star surface, where the redshift is z = (1 — 2GM/Rc 2 )~ 1 — 1. For example, 
T e ff,oo — T e ff/(1 + z) and = F j (1 + z) 2 . As a result, the so-called radiation radius Roo, a 
quantity that can be estimated if F^, T e ff )00 and d are known, is defined to be R^ = R(l + z). 
Roo is a function of the mass and radius of the neutron star, but if redshift information is 
available, perhaps from spectral lines, M and R could be separatel y d etermined. Indeed, 
observation of spectral lines has been reported from IE 1207.4-5209 |48j and EXO 0748-676 
[3], but the identifications of the lines are controversial 0| with redshifts ranging from 
0.12 to 0.35. 

A serious hurdle in the attempt to determine R^ and T e ff tOQ is the fact that neutron 



stars are not blackbodies 
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The star's atmosphere rearranges the spectral distribution 
of emitted radiation. Although models of neutron star atmospheres for a variety of com- 
positions have been constructed, these are mostly limited to non-magnetized atmospheres. 
Pulsars, however, are thought to have magnetic field strengths on the order of 10 12 G or 
greater j^J . The behavior of strongly magnetized hydrogen is relatively simple, but models 
of magnetized heavy element atmospheres are still in a state of infancy 5^] . 

A useful constraint on models is provided by a few cases in which the neutron star 
is sufficiently close to Earth for optical thermal emission to be detected (distinguished by 
green boxes in Fig. 4). These stars have optical fluxes several times less than what a 
blackbody extrapolation from the observed x-rays into the Rayleigh- Jeans optical domain 
would imply. This optical deficit is a natural consequence of the neutron star atmosphere, 
and results in an inferred R^ greater than that deduced from a blackbody. In most cases 
a heavy-element atmosphere adequately fits the global spectral distributions from x-ray 
to optical energies while also yielding neutron star radii in a plausible range. However, 
the observed absence of narrow spectral features, predicted by heavy-element atmosphere 
models, is puzzling The explanation could lie with broadening or elimination of 

spectral features caused by intense magnetic fields or high pressures. 

Radius estimates from isolated neutron stars, while falling into a plausible range, are also 
hampered by distance uncertainties. Pulsar distances can be estimated by dispersion mea- 
sures , but these have uncertainties of 50% or more. In a few cases, such as Geminga fsij , 



HQ and PSR B0656+14 



RX J185635-3754 [5J, m and PSR B0656+14 |£2fl, parallax distances have been obtained, 
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but errors are still large. 

The recent discovery of thermal radiation from quiescent x-ray bursters (involving neu- 
tron stars in binaries) in globular clusters is particularly exciting. At first glance, it seems 
strange that neutron stars in globular clusters, which are on the order of 10 billion years 
old, could be hot enough to emit observable thermal radiation. However, it is believed that 
recent episodes of mass accretion from their companions has been a literal fountain of youth, 
replenishing their reservoir of thermal energy |60|. The measurements of radii from these 
stars might become relatively precise, especially if the distances to the globular clusters in 
which they are found can be refined. Values of in the range of 13 to 16 km have been 
estimated from the quiescent x-ray sources in the globular clusters NGC 5139 and 47 Tuc 



|fil 
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Theoretical cooling curves can be compared to observations if ages for the thermally- 
emitting neutron stars can be estimated (Fig. 4). The best-determined ages are those 
for which dynamical information, such as observed space velocities coupled with a known 
birthplace, is available. Characteristic spin-down ages estimated from pulsar periods P and 
spin-down rates P using t s = P/2P j^| are less reliable. In the cases in which both kinds 
of age estimates are available, they are generally discrepant by factors of 2 to 3. 

Theoretical cooling tracks, for a variety of mass, radius and superfluid properties, are 



relatively narrowly confined as long as enhanced cooling does not occur |43J]. These tracks 
are mostly sensitive to envelope composition. When enhanced cooling is considered, cooling 
tracks fall in a much wider range (Fig. 4). Although most observed stars are consistent with 
the standard cooling scenario, a few cases, espcially PSR J0205+6449 in 3C58 for which 
only upper limits to temperature and luminosity exist may suggest enhanced cooling. 
Uncertainties in estimated temperature and ages have precluded definitive restrictions on 
EOS's or superfluid properties from being made. 



Glitches 



Pulsars provide several sources of information concerning neutron star properties. The 
fastest spinning pulsars yield constraints on neutron star radii. Ages and magnetic field 
strengths can be estimated from P and P measurements. An additionally rich source of data 
are pulsar glitches, the occasional disruption of the otherwise regular pulses Although 
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the origin of glitches is unknown, their magnitudes and stochastic behavior suggests they 



are global phenomena 



64l | . The leading glitch model involves angular momentum transfer 



in the crust from the superfluid to the normal component 33[. Both are spinning, but the 
normal crust is decelerated by the pulsar's magnetic dipole radiation. The superfluid is 
weakly coupled with the normal matter and its rotation rate is not diminished. But when 
the difference in spin rates becomes too large, something breaks and the spin rates are 
brought into closer alignment. The angular momentum observed to be transferred between 
these components, in the case of the Vela pulsar, implies that at least 1.4% of the star's 
moment of inertia resides within the crust leading to the M — R limit in Fig. 2. 

However, observations of long-period (~ 1 year) precesssion in isolated pulsars appear to be 



inconsistent with the crustal glitch model 
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Quasi-Periodic Oscillations 

Quasi-periodic oscillators (QPO's) are accreting neutron stars that display quasi-periodic 
behavior in their x-ray emissions. Generally, their power spectra contain a number of features 
the most prominent of which are twin high frequency peaks near 1 kHz, separated by about 
300 Hz. An early interpretation of these peaks, offered in the sonic point beat-frequency 
model implies a relatively large neutron star mass, M < 2 M [67]. This model holds 
that the higher peak frequency is the orbital frequency of the inner edge of the accretion 
disk and that the separation of the peaks is either once or twice the neutron star's spin 
rate, but fails to account for the observed variations in peak separation as a function of the 
lower peak frequency. Therefore, a variety of other models, most but not all based upon 
rotational phenomenona, are under consideration 671. However, none of these models seems 

n 

to be wholly satisfactory in explaining the observations 67J. 



Future Prospects 

Future observations of binary pulsars and isolated neutron stars hold the promise of 
effective constraints on neutron star maximum masses, radii and internal compositions. 
The importance of the nuclear symmetry energy for neutron stars and supernovae has not 
been overlooked by the nuclear physics community. New accelerator experiments, including 
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high-resolution studies of the neutron skin thickness (which is sensitive to the s ym metry 
energy function S v ) by parity- violating electron scattering on Pb 208 , are planned 16 8l. An- 
ticipated studies of extremely neutron-rich nuclei with rare-isotope accelerators [31( will 
probe conditions intermediate between laboratory nuclei and neutron star matter. Planned 
intermediate energy heavy-ion experiments j3| could establish the in-medium properties of 
pions and kaons that are crucial for delimiting the extent of Bose condensation in dense 
matter. Hyper-nucleus experiments jz(| will shed light on strong interaction couplings of 
strangeness-bearing hyperons likely to occur in dense matter. 

A new generation of neutrino observatories also hold great potential for studies of proto- 
neutron star evolution and neutron star structure. Neutrino observations of supernovae, 
validated by the serendipitous observatons of SN 1987A which yielded about 20 neutrinos, 
should detect thousands of neutrinos from a galactic supernova 
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This could yield 



neutron star binding energies to a few percent accuracy and provide estimates of their 
masses, radii, and interior compositions, as well as details of neutrino opacities in dense 
matter. Neutrino fluxes from proto-neutron stars with and without exotica (hyperons, Bose 
condensates and quarks) have been investigated in jisi I2I | . 

Gravitational radiation is expected from asymmetric spinning compact objects, from 
mergers involving neutron stars and black holes, and from gravitational collapse supernovae 
Depending on the internal viscous forces in rotating neutron stars, gravitational radi- 
ation could drive an instability in r-modes of nonradial pulsations to grow on a time scale 
of tens of seconds 

Q- Mergers Q can be observed to great distances. Detectors due to 
begin operation over the next decade, including LIGO (Laser Interferometer Gravitational- 
Wave Observatory), VIRGO (Italian- French Laser Interferometer Collaboration), GEO600 
(British-German Cooperation for Gravity Wave Experiment), and TAMA (Japanese Inter- 
ferometric Gravitational- Wave Project) could see up to hundreds of mer gers per year 
Binary mergers can yield important information, including the masses |7_3j and mass-to- 
radius ratios of the binary's components and possibly details of their inspiraling orbits jz^. 
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TABLE I: Neutron Star Mass Measurements (la uncertainties) 
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Cen X-3 
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1 - dJ -0.42 


[85, 86] 
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Neutron Star - 
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Neutron Star - Main Sequence Binary 

[22] 



* Could possibly be a black hole, due to lack of pulsations. * Data from 
Reflects binary period- white dwarf mass constraint from jo7|. 
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FIG. 1: The main stages of evolution of a neutron star. Roman numerals indicate various stages 
described in the text. The radius R and central temperatures T c for the neutron star are indicated 
as it evolves in time t. 
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FIG. 2: Mass-radius diagram for neutron stars. Blac k (g reen) curves are for normal matter (SQM) 
equations of state [for definitions of the labels, see [27]]. Regions excluded by general relativity 
(GR), causality and rotation constraints are indicated. Contours of radiation radii are given 
by the orange curves. The dashed line labeled AI/I = 0.014 is a radius limit estimated from Vela 
pulsar glitches j^ . 
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A NEUTRON STAR: SURFACE and INTERIOR 




FIG. 3: The major regions and possible composition inside a normal matter neutron star. The top 
bar illustrates expected geometric transitions from homogenous matter at high densities in the core 
to nuclei at low densities in the crust. Superfluid aspects of the crust and outer core are shown in 
insets. [Figure courtesy D. Page.] 
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FIG. 4: Observational estimates of neutron star temperatures and ages together with theoretical 
cooling simulations for M = 1.4 M . Models (solid and dashed curves) and data with uncertainties 
(boxes) are described in [4^ . The green error boxes indicate sources from which thermal optical 
emissions have been observed in addition to thermal x-rays. Simulations with Fe (H) envelopes are 
displayed by solid (dashed) curves; those including (excluding) the effects of superfluidity are in 
red (blue) . The upper four curves include cooling from modified Urea processes only, the lower two 
curves allow cooling with direct Urea processes and neglect the effects of superfluidity. Models for- 
bidding direct Urea processes are relatively independent of M and superfluid properties. The yellow 
region encompasses cooling curves for models with direct Urea cooling including superfluidity. 
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